Weak Localization Thickness Measurements 
of Si:P Delta-Layers 
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We report on our results for the characterization of Si:P delta-layers grown by low temperature 
molecular beam epitaxy. Our data shows that the effective thickness of a delta-layer can be obtained 
through a weak localization analysis of electrical transport measurements performed in perpendicular 
and parallel magnetic fields. An estimate of the diffusivity of phosphorous in silicon is obtained by 
applying this method to several samples annealed at 850 °C for intervals of zero to 15 minutes. 
With further refinements, this may prove to be the most precise method of measuring delta-layer 
widths developed to date, including that of Secondary Ion Mass Spectrometry analysis. 



The study of delta-layer (5-layer) structures is a very 
active area of research in contemporary semiconductor 
physicsQ, primarily due to the device applications they 
offer, including ultra-high efficiency CCDs, low-energy 
particle detectors^ , and resonant interband tunnel 
diodes (RrTDs)@ . In addition, recent proposals for the 
construction of a silicon-based quantum computer^ , in 
which accurately positioned phosphorous donors function 
as quantum bits (qubits), make the study of Si:P (5-layers 
an obvious first step in the direction of achieving this 
goal 5] . This letter will demonstrate that an analysis 
of weak localization signals in perpendicular and paral- 
lel magnetic fields provides a very precise method for 
measuring the thickness of (5-layers. This quantity is ex- 
tremely important for characterizing the quality of these 
structures, since broadening can occur during growth 
(segregation) as well as in subsequent processing (dif- 
fusion). At present perhaps the most common method 
for determining the thickness of (5-layers is through the 
use of Secondary Ion Mass Spectrometry (SIMS). Un- 
fortunately, SIMS has several drawbacks which limit its 
usefulness in this context. One is that SIMS is a destruc- 
tive technique. More importantly, the spatial resolution 
of SIMS is limited^] to approximately 5 nm. A more sen- 
sitive probe is necessary for applications requiring ultra- 
thin (near monolayer) distributions. Our technique is 
free of these limitations. 

When the two-dimensional electron gas (2DEG)0 as- 
sociated with a (5-doping distribution is cooled to low 
temperatures (T< 4.2K), the resistance of a device fabri- 
cated from such a structure decreases in an applied mag- 
netic field (B). This behavior is characteristic of the phe- 
nomenon of weak localization Q (WL), a quantum effect 
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due to the Aharanov-Bohm interference between oppo- 
sitely directed and closed electron paths in the device. 
At B=0, the amplitudes associated with these paths are 
coherent, interfere constructively, and make the device 
more resistive. For B^O, the interference is destructive, 
and the sample is less resistive. An applied parallel B is 
consequently an extremely sensitive probe of the effective 
thickness of the electron system under study. An early 
WL study on Si MOSFETs (Metal Oxide Semiconduc- 
tor Field Effect Transistors) in parallel fields measured 
the root-mean-square (rms) thickness of the 2DEG with 
sub- monolayer resolution (~ 0.2 nm). If such resolution 
could be obtained in <5-layer studies, it would better SIMS 
capabilities by roughly an order of magnitude. With the 
caveat that the sample must be electrically conducting 
at low temperatures for this technique to work (meaning 
for Si:P a doping density greater than the metal-insulator 
transition (MIT) value of ~ 4x10 18 (cto~ 3 )), we see no 
fundamental reason why monolayer resolution should not 
be achievable. 

A SIMS profile from the wafer used in our experi- 
ments is shown in Fig. 1. This sample was grown by 
low temperature molecular beam epitaxy in a six step 
process. The P <5-layer itself was grown at 320 °C, 
"sandwiched" between Si grown at higher temperatures. 
This (5-distribution has a full width at half maximum 
(FWHM) of approximately 5 nm, close to the resolu- 
tion limit of SIMS. For electrical measurement, four 50 
fxxsx x 50 fjxa van der Pauw devices were fabricated us- 
ing standard lithographic techniques, three of which were 
subjected to an 850 °C anneal for 5, 10 and 15 minutes, 
respectively, to broaden the <5-layer by thermal diffusion 
(henceforth these devices will be referred to by annealing 
time as cook 0, cook 5, cook 10 and cook 15). These 
samples were then mounted in a dilution refrigerator fit- 
ted with a tilting stage which enabled adjustment of the 
angle (9) between the sample plane and B. The base 
temperature of this refrigerator was 130 mK, and from 
temperature-dependent resistivity measurements we esti- 
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FIG. 1: A SIMS profile of an Si:P J-layer sample. The J-layer 
is located approximately 50 nm below the surface, and has a 
FWHM of ~ 5 nm. 



mate that the relevant electronic temperature was ~200 
mK. Resistivity measurements were performed using a 
resistance bridge, four-wire measurement techniques and 
the appropriate Van der Pauw transformation. The sig- 
nals induced in the samples due to a current bias (~ 
1A) were preamplified and recovered using low-frequency 
lock-in detection. We determined the orientation of our 
samples through careful measurements of the Hall coeffi- 
cient, p xy = Bsin(#)/ne (n is the sample carrier density, 
e the electron charge). By taking successive Hall traces 
very close to parallel field we were able to determine 9 
to better than 0.05° (this stray B would contribute only 
~0.2 nm to our thickness estimates determined below). 
In addition, our relatively small sample size minimized 
the influence any B inhomogeneities may have had on 
our measurements. 

For a two-dimensional system the WL signal depends 
on both the magnitude and direction of the applied B. 
When B is perpendicular to the (5-layer (Bj_), the change 
in conductance for a system with weak spin and spin- 
orbit scattering (as is the case for Si:P) is given bv|10| 



Sa(B ± ) = ( 



K 2 



2eB±L 



4eB±L<f, 



(i) 



Here VP is the digamma function, the dephasing 
length (the mean distance over which the wave function 
loses phase coherence), and L the mean free path. In 
parallel fields (B||) the change in conductance becomes 
a logarithmic correction, quadratic in By \H\ and charac- 



FIG. 2: Weak localization data (points) with fits (solid lines) 
to equations (1) and (2) in perpendicular and parallel mag- 
netic fields. The conductance changes much more rapidly 
with applied field in the perpendicular configuration due to a 
larger magnetic flux through the device. 



terized by a single parameter 7: 



M#||) = (^)ln(l + 75f). 



(2) 



The data and theoretical fits for cook are shown in 
Fig. 2, along with the parallel field plots for the annealed 
samples. The difference between the Bj_ and B|| config- 
urations can be understood by the significantly reduced 
flux penetrating the (5-layer when B and the current are 
coplanar. 

To extract an rms 2DEG thickness from the above fit- 
ting parameters, we employ a theoretical model devel- 
oped for MOSFETs [11] in the lowest subband (as far 
as we are aware, a model of weak localization in (5-layers 
does not exist at this time). Of particular concern is how 
the spreading of the donor distribution affects the mag- 
netoresistance. As described in Ref. 12, a surface can 
be characterized by an rms amplitude (here correspond- 
ing to, T, the 5-layer thickness), and a length L c over 
which the fluctuations in the surface height arc corre- 
lated. Different thickness estimates are obtained depend- 
ing on the relative magnitudes of L, L^, and L c , with 
analogs to homogeneous (L C <CL^,) and inhomogeneous 
(L C ^>L^) broadening in nuclear magnetic resonance. For 
(5-layers, the roughness likely arises from fluctuations in 
the positions of individual donors, suggesting a short cor- 
relation length of order the mean donor spacing L c = 
1/i/n (~1 nm). We have no way to independently mea- 
sure L c , so we make this approxiamtion and thus have 
the inequality (see Table I) L c <cL<cL^. This leads to an 
estimate of the (5-layer thickness given by (see Eqn. 55 
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TABLE I: Experimental results for 5-layer samples annealed 
("cooked") at 850 °C for intervals of to 15 minutes. The 
parameters L and L$ were extracted from fits to Eqn. (1), 
7 from fits to Eqn. (2). The mean-square thickness of a 
given layer is proportional to the parallel field parameter 7. 
The dominant error in the thickness estimate comes from the 
mean free path, L. 
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ofRef. 12): 



From this we calculate the diffusion coefficient of Si:P at 
850 °C to be D = 1.6xl(r 14 {cm 2 /sec). Another study 
of Si:P using SIMS measurements reports that for 
this temperature D — 4.16xl0 -16 (cm 2 /sec). This dis- 
crepancy may stem from a number of sources, including 

4000- 
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FIG. 3: A plot of the mean square 5-layer thickness versus 
annealing time at 850 °C. The slope of the line fit to the 
annealed data points yields an estimate of the diffusion coef- 
ficient of Si:P at this temperature: D = l.6xW~ 14 (cm 2 /sec). 
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Using Eqns. (l)-(3) we have calculated the rms thick- 
ness of the S- layer in each of our devices, and shown our 
results in Table I. Our expectation was that the mean 
square thickness of a sample would be directly propor- 
tional to the time it was annealed. Table I clearly estab- 
lishes this relationship, thereby verifying that this tech- 
nique does indeed measure <5Tayer thicknesses. In Fig. 2 
we observe that the parallel field parameter 7 (and there- 
fore the rms thickness) systematically increases with an- 
nealing time. Using these results we can estimate the dif- 
fusivity, D, of Si:P at 850 °C. Assuming an initially Gaus- 
sian profile with thickness To, the mean-square thickness 
after annealing for a time At is expected QjJ to be 
T 2 = T 2 + 2D At. Our data is plotted in Fig. 3. We have 
rejected the cook data point from our linear fit, since it 
is well known that defects such as interstitials and vacan- 
cies enhance diffusion^) . Considering the low tempera- 
ture under which our (5-layer was grown (320 °C), it seems 
likely that a large number of defects were present prior 
to annealing, and that this disorder was subsequently 
" cooked out" (many examples of this sort have been dis- 
cussed in the literature 15]). Fitting a line to the an- 
nealed data points, we obtain a slope of 186(nm 2 /min). 



the thickness model we have employed being inadequate 
(particularly 1/y/n underestimating L c ), different sample 
preparation techniques or perhaps systematic effects as- 
sociated with B (such as imperfect sample alignment or 
field inhomogeneities) . 

Based on our diffusivity results and SIMS data, it 
seems clear that we have systematically overestimated 
the thicknesses of the ^-layers in our samples. We have, 
however, shown a monotonic increase in a devices rms 
(5-layer thickness subsequent to annealing. Therefore we 
believe that by doing more extensive measurements and 
eliminating the various systematic errors in our exper- 
iments (whose contribution to our thickness estimates 
we have not included), this technique will provide the 
most precise method of measuring the thickness of very 
thin ci-layers yet developed. For the future, we plan to 
investigate samples with densities near the MIT limit, 
which should give larger WL signals and perhaps be more 
tractable from a theoretical viewpoint. 
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